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1 The Korteweg— de Vries solitonic spin chain 

i , In a series of papers 0[1O], we described a direct relation between soliton solutions of integrable 

hierarchies and lattice gas systems (e.g., Coulomb gases on two dimensional lattices). The latter 
q models can be reformulated also as some Ising spin systems with a non-local exchange. In 

O . particular, the grand partition functions of specific iV-site Ising chains for some fixed values 

of the temperature were shown to coincide with the tau-functions of iV-soliton solutions of the 
Korteweg-de Vries (KdV) and Kadomtsev-Petviashvili (KP) equations. 

We would like to complete here the consideration of [9] and investigate a critical phenomenon 
£*"»« ', appearing in these models in the zero temperature limit. The iV-soliton solution of the KdV 

equation ut + u xxx — 6uu x = has the form [TJ [11] 

^ i u(x,t) = -2d 2 x log t n (x,t), 

, where tn is the determinant of a N x N matrix M, 



r N = det M, Mij = 5 %J + hB.^)/* 

5-! ' 



k{ H - kj 



0i = hx - kft + ef ] , i,j = 1,2,..., N. 

The parameters k{ describe amplitudes of solitons, 0- /k^ are the zero time phases of solitons, 
and kf are their velocities. The tau- function r^r admits the following Hirota type representa- 
tion 0: 

TN = ^2 eXP ^2 -^Jl'^'j 

;Ui=0,l \l<i<j<N 




*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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where the soliton phase shifts Aij are expressed in terms of the spectral variables ki as 



As remarked in [9] for 6i = 9^°\ this tjv defines the grand partition function of a lattice gas 
model with the chemical potential 6^ and /Xj being the filling factors of the lattice sites by 
molecules. The constants A^ describe interaction energy of the molecules. 

Substituting in (pQ) /Xj = (<7j + l)/2, where crj = ±1 are other discrete variables, one can pass 
from the lattice gases to Ising spin chains [3]: 



t n = e 



1 1 N 



4^ " J ' 2 

i<j j=i 



where 

^=£e-^, (2) 



<Ti=±l 



= E J^o-jcrj- - y^ HjUj. (3) 

l<i<j<N i=l 

Here, Jjj are the exchange constants, ffj is an external magnetic field, T is the temperature, 
and k is the Boltzmann constant: 



N 

A 



4 v ' 2 4 

r^v coincides thus with the partition function of a one-dimensional Ising chain with the specific 
non-local exchange. A similar situation holds for the KP hierarchy and some other partial 
differential or difference nonlinear integrable equations. 

From the thermodynamic point of view, it is interesting to understand the N — > oo behaviour 
of these "solitonic" statistical mechanics models. In general there are infinitely many free 
parameters, and it is difficult to classify qualitatively different cases. An interesting class of 
models is related to the so-called self-similar potentials [12] . These potentials are characterized 



by the (7-periodicity constraints kj + M = qkj and = @j , where q is an arbitrary parameter, 

< q < 1, and M is a positive integer. Exchange constants satisfy in this case the constraints 
Ji+M,j+M = Jij- For M = 1, this translational invariance takes the simplest form J,,- = J(i — j) 
with ki forming one geometric progression, ki = k\q % ~ 1 , q = e~ 2a , where k\ and a > are free 
parameters. More precisely, 

= 2 log | tanh a{i — j) | . 

The KdV coordinate x and time t describe a part of the magnetic field Hi decaying exponen- 
tially fast for i — > oo because q < 1. Only the values of constants 9^ are therefore relevant for 
the leading asymptotics of the partition function in the N — > oo thermodynamic limit. Neglect- 
ing this (x, independence, we come to the constraint Hi + M = Hi, which is just the homogeneity 
condition for M = 1. In principle, it is possible to compute the N — > oo asymptotics for the 
partition function for arbitrary M-periodic magnetic fields, but in [9] only the M = 1,2 cases 
were considered. 

Since < | tanha(i — j)\ < 1, we have Jij = —Aij/A(5 > 0, which corresponds to an antifer- 
romagnetic Ising chain (a similar picture holds for M > 1). Although we have a long distance 
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interaction, its intensity falls off exponentially fast, and the absence of phase transitions in such 
systems at non-zero temperatures is well known [5]. It appears, however, that there exists a spe- 
cial limit leading to a nontrivial critical phenomenon in this model. Indeed, we consider the limit 
a — > + or q — > 1~. The phase shifts A« oc J^jhT then diverge. We can take nevertheless as 
the true exchange constants JJJ n = Jij(q —q) and as the true temperature kT mn = kT(q~ 1 —q). 
For the self-consistency of the temperature definition, we should renormalize the magnetic field 
as well, H oc h/(q~ 1 — q), and assume that h is finite. 

As a result, the interaction energy of any spin "in the bulk" with all others, 



oo 



Ei= ^2 J ij n(J j> 



is finite for q — > 1 (or a — > + ). Indeed, the maximal value of this interaction energy is 



2 



2 

£ max =lim4a V J 0j = lim a V log 

j=— 00,7^0 j=— oo,^0 



l + ql 



4 , (?;?) c 

lim a log ■ 



/3a^o 6 (-<?;?) 



oo 



oo m 

We use the notation (a; 5)00 = fl (1 ~~ a{ ? fe ) an d ( a i> • • • > a m! 9)00 = Tl ( a j! 1 

fc=o j=i 
For the following considerations, we need theta functions [2] 



?2 



73 



nGZ 

(„, 9 ) = ^^+V 2 ) 2 e (2n + 1)„ = ^ e ^_^_^ e -^.^ 

(u,q) = £<fV^ = {q\-qe^\-qe-^^)^ 

n&L 

9,(u,q) = J2(-l) n q n2 e 2nw = {q 2 , qe 2i » , qe~ 2 ™ ; q 2 ) ^ 

where q = e mT , Im(r) > 0, and their modular transformations 

9 1 (v/r,q) = -iy/^iPe^^O^q), 9 2 (u/r,q) = V^Pe iv2 ^ T 9 4 (u, q), 
9 3 (u/r,q) = V^Fe iv2 ^ T 9 3 (u,q), 9 4 {u/r,q) = ^e iv2 ^ T 9 2 (u,q), 

where q = e~ m l T and y '—it is positive for purely imaginary r. Using these formulas, we obtain 
(q;q)l 9[(0,qV 2 ) 9[(0,q^ 2 ) , /8 2 (q; q)^ 



(-q;q)lo 9 2 (0,q 1 / 2 ) (-ir)0 4 (O, Q {-it) (q^-q) 

where q = e 2mT = e~ 2a and q = e~ 2m l T = e~ 2n2 ^ a . As a result, 

"max — 2^ ^ 

and Ei < E mSuX . The limit a — > corresponds thus to an infinitely small and infinitely long-range 
nonlocal interaction model at a low value (zero) effective temperature. 

There are some other interesting limits. For instance, for q — > and finite h, we obtain the 
high temperature nearest neighbor interaction spin chain, JJ,? n oc Si+ij, T ren — > 00. For finite H, 
this limit corresponds to the non-interacting spins. The solitonic interpretation describes thus 



2 ' 

00 
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only a two-dimensional subspace of parameters (T, H, q). For fixed q, the temperature T is also 
fixed, and we can set the "KdV temperature" equal to f3 = 1. 

Using the Wronskian representation for rjy, the leading asymptotics of Zn for N — ► oo was 
determined in [9] for the M = 1 translationally invariant model and a homogeneous magnetic 
field. Namely, — > exp(— iV (3fi), where the free energy per site // has the form 

-Pfj(q, H) = log 2 (g 4 ;g 4 )oocosh/3ff + 1 ^ ^ ^ 2 _ 

(q 2 ;q 2 ^ •* 



(q^^e-^q^ 4sm 2 (v/2) H 6f(u/2,q 2 )' 
The total magnetization of the lattice takes the form: 

N 



m{H) lim 1 V (a.;) = f 1 - - T ^ „ ^ tanh/?ff, (4) 

v ; dH n^oo N 11 V Wo 1 + d(is) cosh 2 PH J w 

where 

g 4 >,g 2 ) 

d W = W( 2T - L 

6>f0, g 2 ) 

The function m(H) grows monotonically with H and reflects qualitative predictions of the 
general theory of ID systems with the fast decaying interactions [5j. However, the limit a — > 
with the renormalized exchange and magnetic field breaks down the corresponding necessary 
conditions, and we obtain a non-trivial critical phenomenon. 

We substitute in (jlj) j3H = ^/(g" 1 — q), h > 0, and take the limit a — > 0. Since 

^ <7 2 ) 0! Wr, <f) (-e 2 -^, -«f e" 2 ^; g 4 ) 2 ^ 



2 (v,q 2 )~ e 2 {v/T,q 2 )~ (e^/r^e^iu/r.^ 
where q = e mr l 2 = e~ 2a , q = e _7n / 2r = q-^ 2 /^> a i we have 



j h {l + e 2lu l T ) 2 {l + e 2i ^~^/ T ) 2 \e h / 2a 



xM = lim cosh — = lim . — —, r— - — 1 

y ' a^o y ' 4a a^o I (1 - e 2w / T ) 2 (l - e 2%( y- li )l T ) 2 j <± 

= lim ( e~ 2iv l T + e 2l( - u -^/ T ) e h l 2a , < i/ < vr. 
Substituting this result in (j3|) and using relation lim tanh(/i/4a) = 1, we obtain 

a— >0 

2 r /2 cfc/ 

m(nj = 1 



K Jo 1 + xO) 
For < i/ < 7r/2, we have 



0, if /l/7T < v < 7r/2, 

oo, if < i/ < h/ir, 



for /i < vr 2 /2, and x(f ) = oo, for h > ir 2 /2. The final result therefore can be represented in the 
form 



m(h) 



2 r 71 / 2 2 

1 - - / dz^ = -zh, if < vr 2 /2, 

1, if \h\ > TT 2 /2. 
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We have an obvious point of non-analyticity of m(h) or of the free energy fi(h) at the cri- 
tical value of the magnetic field /i cr ;t = vr 2 /2, such that the magnetic susceptibility x(h) = 
[3~ l dm(h) / dh has a jump at it (i.e., we have the phase transition of the second order). This is 
a typical phenomenon in the systems with long-range interaction, where the mean field approxi- 
mation gives exact values for the one-point correlation functions (see, e.g., [3]). 



2 The mean field approximation 

In the mean field theory, one considers a few degrees of freedom (usually, just one) of a taken 
system in an effective mean field of the remaining part of the system. This effective or mean 
field depends itself on the analysis of the one-body dynamics. As an example, we consider the 
general spin chain with energy (|3|) and the mean magnetization at the i-th site of the lattice 



-PE 



E a i e ~ 
a\,<72, •■■ 



E e-^ • 

fl,CT2,... 

Instead of calculating the above sums, we stick to the i-th spin and evaluate its contribution to 
the energy as 

Ei(ai) = -aiHi, (5) 

where 

Hi = -J2 J v( a j)+ H i ( 6 ) 

is an effective mean magnetic field at the i-th site created by the external field Hi and the rest 
of the system oc ^ Jij{&o)- 

In the one body problem ([5]), the configuration space consists of two states Oi = ±1, and 
therefore 

^fiHi g — P Hi 

^ = ~m iw = tanh (/^)- 

gpHi _|_ g— pHi 

Substituting the values of effective fields ([6]) in the last equation, we obtain a system of tran- 
scendental equations for mean values of all spins 

(oi) = tanh/3 ^- Jijiaj) + H^j . (7) 

We consider now the translationally invariant system Jij = J(i — j), Hi = H. In such a system, 
all mean values of the spins are the same in the thermodynamic limit, (crj) = (a), and from ([7]), 
we obtain 

(a) =tanh/3 I -J^ J(i-j)(a) +H j = tanh/3(- J (a) + H), J = ^J(j). 



The solution of this equation is an intersection of graphs of two functions: y = x and y = 
tanh/3(— Jx + H). We consider two different cases. 



1. 
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1) A ferromagnet in the zero magnetic field: J < 0, H = 0. Our system of equations 
has only the trivial solution x = y = for (3\J\ < 1 and three solutions x = y = and 
x = y = ±m for j3\ J\ > 1 and some < m < 1. The latter nontrivial solutions describe the 
spontaneous magnetization m at the temperatures smaller than the critical value l//3 CT it = l^l 
(the tanh(— /3Jx)-function becomes steeper at the origin as (3 increases and starts to intersect 
the line x = y in two additional points as its slope exceeds the critical value). 

2) An antiferromagnet, J > 0. The only solution at H = is the trivial solution x = y = 0. 
We consider now the zero-temperature limit (3 — > oo. In this case, the tanh- function transforms 
to the sign- function, and we obtain 

y = x, . . _ j 1, x > 0, 



/ t tt\ sgn (x) — i 
y = sgn(— Jx + H), 6 v ' \ -1, x < 0. 

Solving these equations is rather easy. The function sgn(x) is shifted by H/J from the origin 
along the x-axis. When \H/J\ < 1, it intersects with the line y = x by its vertical part, and the 
magnetization equals to H/J. When \H\ exceeds J, the line y = x intersects with one of the 
horizontal branches of sgn(x), and the magnetization becomes equal to ±1. The magnetization 
is thus a continuous piecewise linear function of H consisting of three parts: two constant 
(a) = ±1 for \H\ > J and the linear piece (a) = H/J connecting them through the origin. In 
our KdV-solitonic model, we denoted (3H = h/(q~ 1 — q) and E m£ix = (q^ 1 — q)J . In the limit 
q — > 1, we have E ma _ x = vr 2 /2/3 and (a) = H/J = 2h/ir 2 , which leads to the exact value of 
the critical magnetic field h cvl t = 7r 2 /2. Such a qualitative behaviour of the system is obvious: 
when the external field exceeds the interaction energy between spins, they all flip in the field 
direction. 

The mean field approximation is known to give exact one-point correlation functions (e.g., 
the magnetization) for systems with the long-range interaction (as our q — > 1 limit). It might 
be non-suitable, however, for the two point correlators (e.g., (fjjCTj)). 

We consider now the M-periodic chain. For general (not necessarily solitonic) Af-periodic 
chain, it is reasonable to introduce multi-index exchange J nm (i — j), where i — j is the distance 
between the cells and 1 < n, m < M are the respective internal cell indices for the n-th and 
m-th sublattices. In our particular solitonic KdV case, we have 



(3J nm (i) = -llog 



2 



J nm (i) = J mn (-i), (8) 



and the energy 



i,j&Z,n^m i,j£Z,ij^j,n i£Z,n 

In the latter sum the magnetic field is also M-periodic, H^ = H^ n \ inhomogeneous only inside 
the cells. The analysis similar to the M = 1 case yields from (|7|) the following system of M 
equations for M unknowns (a^): 



(a (n) ) =tanh/3 ^- ^ </ nrr > (m) > + # (n) ^ 



n, m = 1, . . . , M. 



In these equations, 



jnn = £ jnn^^ jnm = jnm^ n + m. 

i£Z,^0 i&Z 
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In the zero-temperature limit, we have 

(a (n) ) = sgn ( - J nm (<? {m) ) + H {n) 



n, m 



(9) 



where sgn(x) is the sign-function. 

In the solitonic case (El) , the MxM matrix J = J nm 



J' 



' ieZ,i^0 



1 - q l 



1 + 



-A, 



J' 



is symmetric with the constant diagonal: 

k n q k r . 



J 



iez 



(10) 



For M = 2, we have J 12 = J 21 = —B and 

= sgn(,4(a {1) ) + B(a {2) ) + H w ), 
(a^) = sgn(J3<<7«) + A{v&>) + H®). 

If we take the uniform mag netic field = 



H, these equations become symmetric 



in crW and and have the solution (a) = (c^) = (c^) stemming from one equation 

(a) =sgn((A + B)(a) + H). 

For the completely uniform magnetic field, there exists thus a solution when the spins flip 

simultaneously for both sublattices for sufficiently large magnetic fields. 

Simultaneous phase transition exists for all sublattices in the uniform field = • • • = 

H {M) = R ^ when ^(i)^ = . . . = ^(A/)) and all 

spins in all sublattices are aligned simultaneously 
for a sufficiently large H. As seen from ([9]), such a solution exists, if Ylm J nm are e Q ua l) which 
is certainly true for M = 2 because of the permutational symmetry. But it may be not so for 
M > 2. For instance, in the solitonic case (I10p for M = 3, we have 



]Tj 2m 



(X 



E ^ 

ieZ,^0 

E ^ 

ieZ,^0 

E 

ieZ,^0 



1 



1 + g* 



1 



+E l0 ^ 



iez 



l + q l 
1 - g* 



+ E l0 ^ 



l + q l 



iez 



+E 1q ? 



k\( 


t ~ k 2 


kit 


f + h 


k2C 




k 2 i 


f + h 


k 3 ( 


t - h 


k 3 i 


f + h 



+E l0 ^ 

iez 



+ E l0 ^ 



iez 



+E lo £ 

iez 



kiq 1 


-h 




+ k 3 


k 2 q % 


-k 3 


k 2 q l 


+ k 3 


k 3 q l 


-k 2 




+ k 2 



These three sums are certainly different for < q < 1, and in the limit (3 — > oo (which we 
cannot reach within the solitonic interpretation for q < 1), the magnetization would become 
a piecewise linear function of H of a more complicated form than in the M = 1, 2 cases. In our 
model, however, the zero temperature is reached by multiplication of the above sums by q~ l — q 
and taking the limit q — ► 1~ (or a — ► + ). All three sums become then equal yielding the same 
magnetization as in the M = 1 and M = 2 cases. 



3 The BKP solitonic spin chain 

Another Ising chain model solved in [9] appears from the multisoliton solution of the KP equation 
of B type, i.e. the BKP equation [4]. The corresponding partition function has the same form ([2]), 
where the exchange constants are 

1 a _ (oj - aj )(bi - bj)(ai - bj)(bi - Oj) 



@Jij — — ~Aij, e^ 1 



(a,i + a,j)(bi + bj)(ai + bj)(bi + aj) 
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For ai = bi = ki/2, this model coincides with the KdV-inspired model at the twice lower value 
of the temperature obtained after the change (3 — > 2(3. 

The translational invariance of this spin chain, J« = J(i — j), yields 

at = <T\ k = bq l ~\ q = e~ 2a , 

where we normalize oi = 1 and assume that < q < 1 as before. This gives the exchange 

QJ 1 tanh 2 l) 2 /(fr +1) 2 

P lJ 4 8 coth 2 a(i- j) - (6- l) 2 /(6 + 1) 2 ' 

where the parameter 6 is restricted to three regions (because of the 6 — > 1/6 invariance): either 
— 1 < 6 < — q (the ferromagnetic chain, < 0), or q < 6 < 1 or |6| = 1, 6 / —1 (the 
antiferromagnetic chain, Jij > 0). 

In the thermodynamic limit N —> oo, the free energy per site for the homogeneous magnetic 
field Hi = H takes the form [9]: 

2tt 



at /m X i {q,q,bq,q/b;q)oo 1 [ j , io / m 

- /?// iT = - log — — + — / dv log \2p(v) , 

4 {-q,-q,-bq,-q /6;<?)oo 4vr J 



where 



„(,) = cosh 2/Mf + , .<"^ , C r^L^ik + gg^TM 



(-e^,-ge 4!/ ;g)ooV ( & 1 ,qb;q)oo , HJOO 
Taking the derivative with respect to H, we find the magnetization 

m(H)=(l-- [ dv - — — ) tanh2/3ff, (11) 

v ; V * Jo 1 + d O) cosh / v ; 

where 

0i(0/2,g 1/2 ) #2(^ 1/2 ) 

d(l/) " 2 W^j e^v + tfaqV^-e^v-M^qV*) (12) 

with 6 = e*^. For real we have |6| = 1, the choice <fi = i'y, < 7 < 2a, yields q < b < 1, and 
for (/> = -/r + Z7, we have —1 < 6 < — q. The limit 6 — > 1 describes the magnetization for the 
"KdV-spin chain" at the twice lower value of the temperature. A simple test of this expression 
consists in the choice 6 = — 1 corresponding to the non-interacting spins, Jij = 0. In this case 
d{v) = 1, and we obtain m(H) = tanh/?i7 as it should be for the free system. 

We substitute now q = e~ 2a and H = 2h/(q~ 1 — q) in (jlip and consider the limit a — » + . 
The factor 2 in front of h was chosen for coincidence of this model with the the KdV spin chain 
with the effective replacement (3 — > 2(3 (i.e., the twice lower value of the temperature). We apply 
the modular transformation to theta functions in (1121) and obtain 



d(u)=2 6lW2T ^ 1/2) 2 (u/T,q^) 



where r = ia/ir and q = e -27 ™/ T = e~ 2n2 / a . Denoting 

Aqh 



x(v) = lim d(y) cosh 



Q^O 1 
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we therefore obtain 



7T(Tj . 7T(fj , . 7r(2f — 7r) . h 

x(u) = lim 



-e~2^^1 — e~J (1 — e 5 )e Q 

g a n ' 2a M — 6 a I — f3 a 7T * 2a j — g a I 



In the region (j> = i-y, < 7 < 2a, we have 
sin(vr7/2a)(l - e^-^e^ 



7Tf 



X(i/) = lim 



a^o 2sin((7r/2 - v)y/a) 

Since < 7/a < 2, the sin- factors do not influence the asymptotic behaviour, and for < v < 
7r/2, we find xty) = for h/ir < v < tt/2 and x{ v ) = 00 for v < h/ir. As a result, we obtain 
m(h) = 2h/-K 2 for \h\ < 7r 2 /2 and m(h) = 1 for |/i| > tt 2 /2. This is the same picture as for the 
"KdV-chain", as it should be because the limit a — > assumes the limit 7 — > or b — > 1. 

In a similar way, for = ir + 717, < 7 < 2a, and a — > 0, we find b — ► —1, i.e. the trivial 
situation of free spins. The most interesting behaviour appears in the region < (f) < ir, for 
which we find 

2 r /2 eft/ 



7 1 + xH 
where 



xO) = 

As a result, 



0, if /l/(7T - 4>) <u, 
00, if u < h/(n — 4>). 



m W = ~( ^ ' lf h < 9 - » 

7T(7T — 0) Z 

and m(/i) = 1, if > 7r(-7r — <p)/2. The critical value of the magnetic field, for which we have 
the phase transition, 

vr(7r-0) 

"-crit — ± g ' 

depends explicitly on the parameter of the model (ft, and for = 0, we obtain the previous result. 
The Ising spin systems associated with the multisoliton solutions of integrable nonlinear equa- 
tions provide thus the models with phase transitions already in their simplest one-dimensional 
spin chain realizations. 

It is interesting to analyze consequences of the antiferromagnetic nature of the exchange. We 
take for this the KdV-inspired Ising chain and apply different magnetic fields to the odd, H±, 
and even, H2, sites. The corresponding magnetization for the odd sites sublattice was derived 
in [9]: 



m odd (H 1 ,H 2 ) = = lim -y^(°2j-i) 

aHi p^oo p z— / 



. , tanh/3if 2 1 r du 

tanh/lHi ^ / — . 

cosh pHx tt Jo ^py. - tanh/3^1 tanh 0H 2 
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The magnetization for the even sites sublattice is obtained after permuting Hi and H2 in this 
expression. Obviously, if we take the alternating magnetic field Hi = —H2 = H, then the total 
magnetization is equal to zero, though the sublattice magnetizations remain non-trivial: 



However, after substituting f3H = hj {q~ l — q) and taking the limit a — ► 0, we see that our zero 
temperature critical phenomenon disappears: m dd(h) = 1 for h > and m Q dd(h) = — 1 for 
h < 0, similar to the free spins system. 

4 Conclusion 

As shown in |10j . soliton solutions of integrable hierarchies with the complex values of spectral 
variables are connected to the intrinsic Coulomb gases on two dimensional lattices with some 
nontrivial dielectric or conductor boundaries. In this picture, the Coulomb interaction energy 
between two charges and their effective images created by the boundary conditions plays the 
role of the soliton phase shifts, the coordinates of charges coincide with the spectral parameters 
of solitons, and the external electrostatic field is expanded in some series with the coefficients 
playing the role of integrable hierarchy times. This transparent relation serves as a clue for 
building new Coulomb lattice gas models exactly solvable at some fixed temperatures. The 
latter temperatures are also related to the random matrix models [H EJ [TO], but we do not 
discuss here applications of the described phase transition within these interpretations. 

Our phase transitions are of a rather simple nature. In the lattice gas language, the transition 
in the M = 1 periodic case describes the situation when the lattice is filled to its limit, i.e. the 
number of particles equals to the number of sites, and no more particles can be added to the 
system. In the M-periodic case, such transitions may happen separately, when each sublattice is 
filled completely one by one, or simultaneously, for an appropriate choice of parameters. Their 
qualitative features can be found from the mean field theory. There are several interesting 
questions which would be interesting to analyze in the future, like influences of the hierarchy 
times on the thermodynamical quantities, understanding of our systems beyond the "solitonic" 
temperature values, investigation of the higher correlation functions, and so on. 
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/ vaguely remember a tall man with glasses vigorously explaining to me something during my 
poster presentation at the IGTMP colloquium in Moscow in 1990. Probably that was Vadim - 
I never asked him about that later on. We got acquainted at the first SIDE meeting near Montreal 
in 1994 and had sufficiently long discussions during his few days visit to CRM after that confe- 
rence. I remember telling him that by the work on separation of variables [8], which impressed me 
much, he closed to me that field, and it is necessary to think about other directions of research. 
We became closer during Vadim's stay at the CRM in 1994-1995. It was very nice time from 
many points of view. I visited his house and have known his family during one of the parties 
he was gathering. It appeared that he likes Russian "bards" singing, which I was bond to as 
well. Once he even sent to me a web-link to some new mp3-recordings coming from his native 
Saint-Petersburg. 
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In June 1999, he chaired my talk at the Hong Kong meeting on special functions, where 
I reported results of the paper [TH] ■ He was interested by that much, and we discussed possible 
intersections with integrable systems. I saw Vadim last time at the conference in Edinburgh in 
September 2003, where he was the main organizer. During the preparation of the corresponding 
proceedings, I actively communicated with him, and it was clear that he is extremely busy by all 
kinds of obligations. In May 2005, I suggested to him to form a team in order to try to get an 
INTAS grant. His first reaction was positive, but after he has known the rules and procedures, 
he rejected this idea by saying that there is too much bureaucracy and he has too many other 
commitments for the next few months. I totally agreed with his critics and accepted his excuse. 
Later on he listed to me a number of other possibilities to get research funding from the UK 
sources which sounded quite reasonable. We exchanged by about ten e-mails with him over the 
May-October 2005 period, and it was devastating to know that he has passed away. We have 
another deeply regrettable loss in the FSU scientific community, which was possible, probably, to 
prevent in other circumstances. 

V.P. Spiridonov 

I have met Vadim first when I have been pursuing my PhD studies at the Centre de Recherches 
Mathematiques in Montreal. At that time, Vadim was a postdoctoral fellow there. I remember 
him as an open-heart person, frank and honest, very enthusiastic and completely devoted to the 
problems he did and had in mind. Always full of energy, he showed the keenest interest for many 
questions of science and life. To meet Vadim was very interesting to me. 

I.M. Louts enko 
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